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1. Introduction 

Transformation Optics (TO) [H [2] has attracted great interest in the past few years 
due to its potential abihty to control the flow of light. It establishes a procedure 
that provides us with expressions for the electromagnetic (EM) material parameters, 
dielectric permittivity and magnetic permeability, that need to be implemented in 
order to obtain a medium where EM waves propagate along desired trajectories. It 
is based on the mathematical equivalence between coordinate transformations and 
EM material parameters. As Maxwell's equations are form-invariant under coordinate 
transformations, a distorted geometry in an empty space can be interpreted as a 
medium in a Cartesian flat space [3j. A wide variety of optical devices has been 
designed by means of TO, such as invisibility cloaks [H [4], beam shifters and dividers 
[5], beam bends ^ ^7\^ waveguide bends [8j, fleld concentrators [9], imaging devices 
[ini [H] or polarization splitters and rotators [12j. The permittivity and permeability 
tensors required by the TO technique are usually implemented using metamaterials, 
artificial materials that allow to eflFectively mimic EM responses by repeating structured 
unit cells at a subwavelength scale [131 El- However, these parameters are generally 
anisotropic and inhomogeneous, so metamaterials usually rely on resonating structures 
to tailor EM parameters, a feature that limits the wavelength range of operation and 
causes losses. Conformal and quasiconformal \Xhi coordinate transformations have 
been proposed to overcome this problem, since they reduce anisotropy to a minimum 
yielding the possibility of fabricating transformation media out of isotropic dielectrics. 
Consequently, devices based on conformal transformations can be inherently broadband 
at optical frequencies. While conformal mappings are usually obtained analytically [16j, 
quasiconformal mappings rely on numerical grid generation techniques. Recently, a 
method based on solving Laplace's equation with the appropriate boundary conditions 
has been proposed [d HE] to generate quasiconformal mappings for different optical 
elements such as cloaks or arbitrary waveguides. 

It has been recently shown that TO recipes can be applied to efficiently mould 
the ffow of surface plasmon polaritons (SPPs) flHl [20l EU [22l [23]. SPPs are highly 
confined waves that propagate along a metal-dielectric interface. They are evanescent 
in the direction perpendicular to the interface and their field decays exponentially in the 
metal and dielectric sides. SPPs arise from the resonant interaction of light with the free 
electrons of the metal and result in a collective oscillation [211 ES] . SPPs are considered 
as good candidates for photonic circuitry as they provide large field enhancements at 
subwavelength scales [26l [27j. On the other hand, within the field of Metamaterials, 
plasmonic structures have already been developed to achieve cloaking. They include, 
for instance, the cloaking from an SPP of an object surrounded by a circular array 
of point scatterers [28], the experimental verification of a plasmonic shell capable of 
cloaking subwalength objects [^ or a two-dimensional cloak for SPPs consisting of 
layers supporting plasmonic modes with group velocity of opposite signs [30]. However, 
it is interesting to note that using the TO tools to manipulate SPPs is a very general 
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technique that can yield different optical devices and new functionalities. Furthermore, 
as it has been shown in previous works [191 ^ o^ily the dielectric material next to the 
metal needs to be structured, while modifying the metal properties is not required in 
most cases. 

In this manuscript we consider the main features of TO as applied to SPPs and 
exemplify the theory by means of different TO devices for SPPs. We verify the results 
using the Finite Element Method solver COMSOL MULTIPHYSICS. First, we introduce 
the general methodology of TO in section [2] and apply it to the case of SPPs. In 
particular, we study the functionality of four different potential plasmonic devices 
(a cylindrical cloak, a parallel shifter, a right-angle bend and a lens) at a range of 
wavelengths through the optical range of the spectrum. Second, we consider realistic 
models based on conformal and quasiconformal mappings in section |3| We show that 
for these cases only isotropic refractive indexes need to be implemented and we apply it 
to the parallel shifter, the right-angle bend, the lens and a ground-plane cloak for SPPs. 
Finally, section [4] summarizes the main contributions of our work. 



2. Transformation Optics for Plasmonics 

The TO framework provides tools for the design of EM media with desired 
functionalities. It is based on the fact that Maxwell's equations are left invariant if 
a mapping is performed from an empty space with a distorted geometry (virtual space) 
to a medium characterized by e and jit placed in a flat space (physical space). A two- 
dimensional (2D) transformation = x^{x^') between the virtual system, x^' = {x',^'}, 
and the physical system, = {x, 7/}, is described by means of its Jacobian matrix: 

A -(IS.) (1) 



^dx'^\ 

where the index i = 1,2. As Maxwell's equations are form-invariant, the empty-space 
equations in physical space can be interpreted as macroscopic equations for a medium 
placed in virtual space and characterized by the following EM parameters: 

AA^ . AA^ . . 

^To = ^0 ^; Mto = Mo ^ (2) 

det A det A 

where cq and /llq are the background EM parameters of the physical space. On the other 

hand, we can relate the metric tensor G = {gij) of the coordinate transformation to the 

inverse Jacobian matrix: 

G=$-|^=fl-f A- (3) 



dx^ dx^ 

This allows us to re-express e^o and jiTO in terms of the inverse metric tensor, 
= AA^ ^ of the transformation: 

^TO = e^^/gG~^] jiTo = l^oy/gG~^ (4) 

where g is the determinant of the metric, g = l/(detyl)^. In general, the metric tensor 
of any 2D coordinate transformation is a 2 x 2 symmetric matrix. Then, as it can 
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be readily seen from Q, e^o and jiTo are generally anisotropic. Moreover, they are 
inhomogeneous and vary throughout space according to the coordinate transformation. 

Since TO relies on Maxwell's equations, it applies to all kinds of EM waves and 
in particular to SPPs. However, we need to bear in mind the fact that SPPs are 
surface waves propagating along a metal-dielectric interface with a field distribution 
that is evanescent in the direction perpendicular to the interface. Therefore, an infinite 
number of 2D coordinate transformations should be performed, in principle, in planes 
parallel to the metal-dielectric interface in order to operate over the whole plasmonic 
field. As the SPP field extends into the dielectric and the metal, these 2D coordinate 
transformations should be implemented at both sides of the interface. Due to the fact 
that the SPP decay length in the metal (skin depth) is of a few tens of nanometers 
at optical wavelengths, a manipulation of the metal properties would be needed at a 
nanometer scale. However, it has been shown [19] that a more feasible approach that 
consists of of modifying the EM properties of the dielectric side only indeed gives very 
accurate results, while control over the metal is not actually required. To illustrate the 
direct application of TO to SPPs and the impact of manipulating only the dielectric side, 
we have studied the dependence of the efficiency on the wavelength for four different 
TO devices for SPPs: a cylindrical cloak, a parallel shifter, a right-angle bend and a 
lens. 



(a) (b) 




Figure 1. Transformation Optics is used to cloak a cylinder to an SPP. (a) An SPP is 
propagating along a gold- vacuum interface and scatters on a metallic cylinder of radius 
a. The plasmonic field intensity, shown schematically, decays exponentially away from 
the interface, (b) Scattering losses are suppressed by placing a cloak of outer radius h 
around the cylinder. The optical parameters needed to cloak the cylinder are calculated 
using TO and implemented in a cylindrical shell a < r <h oi height hd in the dielectric 
side {metamaterial 1) and in a shell of height hm in the metal {metamaterial 2). As 
the decay length of the SPP is larger in the dielectric than in the metal and most part 
of the SPP energy is carried in the dielectric, hm is much smaller than hd. 



2.1. Cylindrical cloak 

As a first example of TO applied to plasmonics, let us consider a three-dimensional (3D) 
cylindrical cloak for SPPs travelling along the interface between a metal characterized 
by permittivity and vacuum [19j. Figure [l](a) shows an SPP that scatters with a 
metallic cylinder of radius a placed on the surface. A cloak placed around the cylinder. 
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see figure [T](b) , renders the object invisible by suppressing any scattering from it. A 
2D invisibility cloak was the first proposal of TO applications [U H] as well as the first 
experimental realization |31j- The cloak is designed by taking a physical space that has 
a hole of radius a surrounded by a compressed region of radius h and then making a 
transformation to a Cartesian undistorted space. A wave travelling in physical space 
will follow trajectories that surround the hole so that it appears to have travelled in 
straight lines as it would in the virtual - undistorted - space. In order to design a 3D 
invisibility cloak for SPPs we take a plane parallel to the metal surface and make the 
2D radial transformation r{r') = b/{b — a)r' + a in the shell a < r < b (see the geometry 
of the cloak in the inset of figure [2]( a)). In Cartesian coordinates, the transformation 
reads: 

X rir')— ; y r(r')— ; z = (5) 

According to the TO equation Q, the EM parameters of the transformation medium 
are the following: 

/ e,^ + e4 \ 

eTo/eo = Ato/mo = {er - ee)^ e^^ + e^^ (6) 

V 6, y 

where Cr = fir = {r — a)/r^ eo = fio = r / [r — a) and Cz = jJ^z = [b/ {b — a)]^ {r — a) / r . These 
parameters are highly anisotropic and approach singular values at the inner radius of 
the cloak {exx-> ^yy +oc, e^^^ ^ — oc, ^ 0). If we want to operate on an SPP, we 
need to cloak the field in the dielectric and metal sides as illustrated in figure [l](b) . In 
the dielectric side, we implement the parameters, e = cto and /x = /xto, in a cylindrical 
shell of radii a and b with a height, hd^ larger than the vacuum decay length of the SPP. 
On the other hand, in the metal side the EM tensors need to be modified, e = cto ' 
and ji = jiTO) in a cylindrical shell of the same radii but with a height hm larger than 
the skin depth. If we cloaked its entire field, an incident SPP would emerge at the other 
side as if it had travelled through free space. 

Figure [2] shows the results obtained from 3D simulations. To study the effectiveness 
of the plasmonic device we have considered the transmittance of SPPs through the cloak 
and its scattering properties. We have measured the transmittance from one side of the 
cloak to the other as the ratio of the incident to the transmitted power flow. In order to 
isolate the effect of the TO medium, losses have been neglected throughout this work 
by setting the imaginary part of the metal dielectric constant to zero. Since losses 
are absent, 100% transmittance means that the whole plasmonic field is guided around 
the cylinder without any back-scattering. In the simulations, an SPP at a wavelength 
between 600 and 1500 nm and travelling through a gold- vacuum interface, impinges onto 
a cloak of size a = b/2 = 5/6A. The transmittance spectrum is plotted in figure [2](a) 
for two cases: with and without introducing a metamaterial in the metal side (blue and 
red lines respectively). The transmittance of the cloak in the ideal situation (blue line) 
where TO parameters are implemented in the metal along with the dielectric material is 
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higher than 99% over the whole wavelength range. On the other hand, the transmittance 
when the metal is leaved untouched (red line) increases from 97% at 600 nm to 99% at 
1500 nm. As the wavelength increases, the blue and red curves lie closer. The reason 
for this is that while the skin depth is approximately constant over the whole spectrum 
(around 20-30 nm), the vacuum decay length increases with the wavelength, varying 
from 250 nm for an SPP at 600 nm to 4 /xm at 1500 nm. Thus, for shorter wavelengths, 
the part of the SPP energy carried within the metal is larger. However, although the 
effectiveness of the cloak is affected by the lack of cloaking the field that resides in the 
metal, it is very good even at short wavelengths. The performance of the cloak for 
an SPP at A = 900 nm is illustrated in figure [2](b), which shows how the SPP field is 
smoothly guided around the metallic cylinder. 

Finally, a study of the scattering properties is needed for a complete analysis of 
the invisibility cloak. By inspection of the scattered field distribution we have seen 
that, when the TO recipe is not implemented in the metallic side, around 20% of the 
field amplitude is scattered in the forward direction. As a result, the wavefront of 
the SPP when passing through the simplified cloak is slightly modified. In accordance 
to [32j, phase distortions occur for the simplified cloak and perfect cloaking is not 
achieved. Nevertheless, as seen in figure [2| most of the plasmonic field is guided around 
the cylinder. On the other hand, we have checked the scattering properties of all the 
devices shown in this work. We have found that these eflFects are only relevant in the 
case of the cylindrical cloak. For the sake of brevity, we mention this fact here and 
not when discussing each device. Therefore, in the rest of the paper, we present the 
transmittance and the images of the EM field as criteria for quantifying the performance 
of the SPP-devices. 



2.2. SPP parallel shifter 

Another potential plasmonic device is the parallel shifter , which acts on an incident 
beam by translating it perpendicularly to its propagation direction. To design a 2D 
shifter, we take a rectangular region sized d x / (see inset panel of figure [3]( a)) and map 
a Cartesian grid to a grid that is tilted at an angle (j) [5] : 

X = ; y = + a{x^ + d) ] z = z' (7) 

where a determines the tilt angle a = tan (j). The map of the transformation is illustrated 
in figure |3](b) . In contrast to the transformation leading to the cylindrical cloak ^ , 
the transformation ([7| is discontinuous at the frontiers of physical space, a feature 
that allows the transfer of the modification of the fields inside the TO medium to the 
wave that exits it. Because the transformed boundary is a combination of rotation and 
translations of the original boundary, the discontinuity does not cause refiections at the 
limits of the TO medium [16j. Following the TO equations ^ the EM parameters of 
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wavelength (nm) 

Figure 2. Simulation results of a 3D cylindrical cloak for an SPP. (a) Transmittance 
(ratio of the incident power flow to the transmitted power flow) through the cloak for 
SPPs at different wavelengths. The blue and red lines show the case with and without 
placing a transformation medium in the metal side respectively. The in-plane size of 
the cloak is scaled with the wavelength a = 6/2 = 5/6A. The height of the cloak in the 
dielectric side, hd^ is much larger than the vacuum decay length at a given wavelength 
and is chosen for convergence. The height of the cloak in the metal side, hm^ is 50 
nm high at all wavelengths since the skin depth of the SPP is very similar over the 
wavelength range and amounts to 20 — 30 nm. The inset panel shows the geometry 
layout from a top view, (b) z-component of the electric field of an SPP at 900 nm 
travelling through a cloak without applying TO in the metal side. 
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the transformation medium are obtained: 

a \ 

(8) 

A 3D SPP parallel shifter should, in principle, be composed of a slab in the dielectric 
side sized dxl x hd with permitivitty cto and permeability fiTO along with a slab in the 
metal film of dimensions d x I x and parameters e = cto ' and jl = firo- Again, 
the height of the slabs at each side must be larger than the SPP decay length in that 
medium. In this case, the TO medium parameters are anisotropic but homogeneous 
and take non-singular values. 

Figure [3](a) shows the transmittance spectrum resulting from 3D simulations. An 
SPP with a gaussian profile of width A = / — 2ad propagates in the x direction and 
impinges on the transformation medium. Thus, it experiences a shift of = 45° (we 
set a = 1) perpendicular to its propagation, in the y direction, and exits the shifter 
propagating along a direction parallel to the original incident path. The in-plane size 
of the device is given by / = 20/3A and d = 1/4. The blue and red spectra in figure |3](a) 
show the transmittance through the beam shifter with and without introducing a 
transformation medium in the metal side, respectively. In the ideal situation (blue 
line), 100% of the SPP powerflow is shifted over the whole spectra, whereas in the 
approximation (red line) the transmission is lowered due to the fact that the fraction 
of the SPP energy that is carried inside the metal is not shifted. At 600 nm the 
transmission is lowered to 95.5% while at 800 nm amounts to 98.5%. On the other 
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hand, for wavelengths larger than 900 nm the red and blue lines differ by less than 1% 
because most of the SPP energy is carried outside the metal. Therefore, manipulation 
of the metal is not required since the performance of the SPP shifter is very good when 
the TO parameters are introduced in the dielectric side in the whole optical range. The 
effect of the shifter on an SPP can be seen in figure [3](c), where the field distribution of 
an SPP at A = 900 nm undergoing a 45° shift is shown. As a consequence of the shift, 
the outgoing propagation direction is 1.5 /xm from the incident path. Note that the SPP 
wave fronts nicely follow the vertical lines of the transformation map shown in (b). 



2.3. SPP right-angle bend 

The next plasmonic element we consider here is a right-angle bend, a device that rotates 
the propagation direction of an incident SPP by | (see the inset of figure |4]( a)). We 
start by describing a 2D beam bend, which is designed by transforming a a squared 
region of side 6 in a Cartesian grid into a polar grid of squared cross-section b x b [7J. 
In order to bend a beam that propagates in the x direction by |, so that it exits the 
device propagating in the y direction, the appropriate transformation is the following: 

r = y'; (j)=^ib-x'); z = z' (9) 

with r = {x'^ + y'^y/'^ and (f) = arctan(^/x). Again, this transformation is not continuous 
at the boundaries of the domain. The EM parameters retrieved from TO equations Q 




Figure 3. 3D simulation results of a parallel shifter. A gaussian SPP of width A 
impinges on the device and experiences a shift of = 45°. (a) Transmittance curve 
of the shifter for SPPs at different wavelengths. The blue and red lines show the case 
with and without placing a transformation medium in the metal side, respectively. The 
TO parameters are implemented in a slab sized d x I x hd in the dielectric side and, 
only for the blue line, in a slab d x I x hm within the metal (see the inset panel for 
geometry parameters). The transmittance through the shift is measured as the ratio 
of the power flow that impinges the entrance face within —A/2 < y < A/2 to the 
power flow that exits within < < A. (b) Coordinate transformation given by 
(c) Top view of the magnetic field component of an SPP at 900 nm travelling through 
a shifter without applying TO in the metal side. Note that the two lines at 45° are 
simply a guide to the eye and do not represent a physical boundary. 
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are expressed in Cartesian coordinates according to ^ with: 

2b nr 2b 

6^ = — ; eo = — ] e, = — (10) 
nr lb nr 

These tensors are singular at r = ( e^, ^ +oc and cq ^ 0). In order to avoid the 

divergence at r = we hmit the extreme values of the parameters by limiting the radius 

to vary from an inner value r = a to the outer value r = b. Then the width of the bend 

is A = 6 — a and the curvature radius is p = {a + b)/2. According to our methodology, 

to design an SPP right-angle bend we need to perform the coordinate transformation 

in the cylindrical shell comprised between r = a, r = 6, = and = |. Ideally, we 

would need to implement the parameters cto and ^to given by ^ and (10) within the 



dielectric material up to a height of hd^ while the metal optical properties should be 
modified following the expressions e = cto ' and fi = fiTo in a region of height /i^. 

Figure [4](a) shows the transmittance curve of an SPP bend for two different 
situations: with and without implementing the TO parameters in the metal side (blue 
and red lines, respectively). The transmittance through the bend is calculated from 
simulations where SPPs at wavelengths ranging from 600 to 1500 nm are launched 
into the entrance of a bend of curvature radius p = |A. The field distribution results 
from a simulation of an SPP at A = 900 nm propagating in the x direction and 
travelling through a right angle bend of radius p = 1.5 /xm is shown in figure [4](b). 
The transmittance is obtained from measurements of the power fiow at the entrance 
and exit faces of the bend. The blue line in the transmittance curve corresponds to the 
ideal case where the TO parameters are implemented both in the metal and dielectric 
sides. The transmission is over 99% for all wavelengths. On the other hand, the red line 
follows the transmittance through the bend when only the dielectric side is modified. 
It is clearly seen that this is indeed a very good approximation, as the transmission is 
always higher than 95%, approaching the same results as the ideal simulations at large 
wavelengths. For wavelengths larger than 900 nm, the blue and red curves differ by 
only 2% or less, meaning that by not structuring the metal side less than 2% of the 
SPP powerfiow will be lost. On the other hand, for short wavelengths (600-800 nm), the 
transmittance of the approximate situation is around 4% lower than the transmittance 
of the ideal blue curve, a value that corresponds to the fraction of the SPP energy that 
is contained in the metal. Therefore, we have shown that in the approximate but more 
feasible situation of structuring only the dielectric side, the performance of the SPP 
bend is very good over the entire wavelength range considered. However, in some cases 
the SPP wave fronts might suffer some distortions even though the transmittance is still 
very high. 

24. SPP lens 

Finally, we consider a plasmonic lens capable of focusing an incident gaussian SPP of a 
finite width. As described in [7J, transforming a semi-circular region of radius a into a 
rectangular region sized ax2a leads to a 2D rectangular lens whose focal point is located 
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Figure 4. 3D simulation results of a SPP right-angle bend, (a) Transmission through 
the bend for SPPs of gaussian profiles (width A) at different wavelengths. The blue and 
red lines show the case with and without placing a transformation medium in the metal 
side, respectively. The bend curvature radius is p = 5/2 A and the width is A = 10/3A. 
The inset panel shows the geometry layout from a top view, (b) z-component of the 
electric field of an SPP at 900 nm travelling through the bend without applying TO 
in the metal side. 



at the center of its exit face (see an schematic layout in the inset panel of figure [5](a)). 
Assuming an incident wave from the left, we place the lens in the x — y plane in the 
region defined by — a < x < and —a<y<a. For such a geometry, the appropriate 
transformation is found to be: 

ax 



X 



y = y;z 
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The TO equations Q yield the following EM tensors in cartesian coordinates: 



I x'^y'^a~^{a'^—y'^)~^-\-a 
xy 



xy 



V 











(12) 



These parameters are inhomogeneous and anisotropic and take singular values at the 
boundaries of the device. As in the previous examples, the lens for SPPs consists of a 
rectangular slab sized a x 2a x hd and placed on top of the metal film, with parameters 
cto and jj^To- In principle, a rectangular slab of dimensions ax2ax hm and parameters 
^TO ' and fiTO should be implemented in the metal side. 

Figure [5] shows the results obtained from 3D simulations of a gaussian SPP of 
width A travelling through the lens. As a result, the SPP is focused at the exit side 
of the device and continues to propagate in cylindrical wavefronts. Panel a shows the 
transmittance through the device as a function of the SPP wavelength for a quantitative 
analysis of the lens effectiveness in two situations: when the TO parameters are or 
are not implemented in the metal side (blue and red curves, respectively). In the 
ideal case where transformation media are placed in the metal and dielectric sides, the 
transmittance is close to 100% at all wavelengths. On the other hand, the red line 
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Figure 5. Simulation results of a 3D SPP lens. An SPP exhibiting a gaussian profile 
of width A and propagating in the x direction impinges on the lens, where it is focused 
to a width of the order of the wavelength, (a) Transmittance spectra through the 
lens for two cases: with and without placing a transformation medium in the metal 
side (blue and red lines respectively). The transmittance is measured as the ratio of 
the power flow incident on the left to the power flow that exits within a wavelength 
width. The in-plane size of the lens is a x 2a, with a = 2.5A, while the width of the 
incident beam is A = 3. 75 A (see inset panel), (b) 3D view of the y-component of the 
magnetic field of an SPP at 900 nm travelling through the lens. In this simulation, 
the transformation medium was placed in the dielectric side only. Note that the lines 
parallel at the exit face separated by 900 nm constitute an auxiliary element to measure 
the transmittance and do not represent a physical boundary. 



shows that the approximation consisting of not modifying the metal side is very good 
for wavelengths above 900 nm as in that wavelength range the transmittance lies above 
98%. For wavelengths below 900 nm the transmittance is lowered up to an 8% compared 
to the ideal case because for shorter wavelengths there is a higher fraction of the SPP 
energy contained in the metal. However, the field patterns at these wavelengths (not 
shown here) reveal that the wave propagation and focusing is still very good. Panel b 
shows a plot of the field for an SPP at 900 nm travelling through the lens. While the 
SPP travels a distance of a = 2250 nm, its wavefronts are smoothly focused on a spot 
of the order of its wavelength. 

3. Conformal and quasiconformal mapping for realistic models 

The coordinate transformations we have considered so far present some technical 
problems when it comes to practical realizations as they generally yield highly 
anisotropic e and fi. However, if the we choose an appropriate coordinate transformation 
by imposing some mathematical requirements, isotropic EM parameters may be 
obtained. This means that only isotropic dielectrics would be necessary to construct 
devices, reducing losses and opening up the possibility of broad-band performances. In 
particular, conformal and quasiconformal transformations give rise to isotropic material 
parameters, yielding easier implementations of transformation media based on isotropic 
dielectrics. Conformal mappings were first investigated in [4], where an invisibility 
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device was designed by applying a transformation in the complex plane to Helmholtz's 
equation. This concept was also applied to design other isotropic devices, such as 
arbitrary waveguide bends within the limit of geometrical optics [33[ [34] or directional 
antennas, flat lenses and beam bends [16j. Regarding plasmonic systems, conformal 
and quasiconfomal transformations have been employed to design of light-harvesting 
[35] devices and ground-plane cloaks [20j. 

In the following we will show how to derive the isotropic EM parameters that 
correspond to a conformal transformation within the TO framework. Let us start by 
assuming that the virtual and physical domains are contained in the complex planes 
z = x' -\- iy' and w = x + iy^ respectively. As we are dealing with SPPs, 2D 
waves extending into the third dimension (the direction perpendicular to the metal- 
dielectric interface), we can restrict our study to 2D transformations. A 2D coordinate 
transformation from virtual to physical coordinates can then be described by means of 
a complex function: 

w = f{z, z) = x{x', y') + iy{x', y') (13) 

This transformation is said to be conformal if the function f{z) depends on z but not on 
its conjugate z and is analytical, i.e. f'{z) can be deflned. Then, the Cauchy-Riemann 
equations are satisfied {dx/dx' = dy/dy' and dx/dy' = —dy/dx'). Consequently, the 
metric tensor of a conformal transformation reduces to a diagonal matrix: 

= . I,., where = (^n)"^ = ^ + ^ (^4) 

This means that the grid described by the transformation map w = f{z) is orthogonal, 
i.e., the angles between coordinate lines are preserved after the mapping. According to 
the TO rule Q, we calculate the EM parameters of the transformation medium that 
corresponds to a 2D conformal transformation: 

/ 1 \ 

^To/^o = P^To/l^o = VgG'^ = 10 (15) 

Now let us consider the case of an EM wave with the electric field polarized along the z- 
axis and the magnetic field in the x-y plane. The propagation of this wave is determined 
by the electric permittivity along the z-axis and the magnetic susceptibility tensor in 
the x-y plane. According to the equation above, these are the following: 

= ^0^9 ; A = Mol2x2 (16) 

This set of parameters is equivalent to an isotropic refractive index: 

n = y/e^i = n^g^'^ (17) 

where no = s/e^jl^ is the background refractive index. If on the other hand, the EM 
wave is polarized with the magnetic field along the z-axis, then e = eg, /x^ = Moy^ and 
n = ^/efiz = '^05'^^^- Notice that, as shown in the Appendix, our approach leads to the 
same result as the Optical Conformal Mapping technique based on Helmholtz's equation 
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[Ij. Therefore, for problems where the two polarizations can be treated separately, a 2D 
conformal transformation leads to an isotropic transformation medium, with a refractive 
index defined by ( 17). On the other hand, SPPs need to be considered separately as their 
dispersion relation depends on the electric permittivity and the magnetic permeability 
rather than on the refractive index. However, we will show in the following that the 
propagation of an SPP in a medium with the exact diagonal anisotropic parameters (15) 
is very similar to that in a medium with the approximate isotropic refractive index (17) 
for a given conformal transformation. 

Thus, conformal transformations lead to transformation media described by 
isotropic refractive indexes. However, an analytical conformal transformation is not 
known for every desired functionality. In those cases, quasiconformal maps are useful 
for designing TO media based on isotropic materials. They constitute a subset of 
coordinate transformations, w = f{z^z)^ that satisfy less requirements than conformal 
transformations. In particular, they can depend on z and z. In contrast to conformal 
transformations, quasiconformal transformations map squares in the original space 
to rectangles of constant aspect ratio in the transformed space. Since they are 
generated numerically through the minimization of, for instance, the Modified-Liao or 
Winslow functionals [36l|37j, they are able to treat arbitrarily complex boundaries. A 
quasiconformal map is the optimal map of minimum anisotropy for a given domain, a 
fact that was can be used to design isotropic transformation media. In reference [15], a 
quasiconformal map for a ground-plane cloak was obtained by minimizing the Modified- 
Liao functional with sliding boundary conditions. Applying this idea, several invisibility 
ground-plane cloaks operating at optical [38l [39] and microwave jlQ] frequencies were 
fabricated out of isotropic dielectric materials and a non-resonant metallic metamaterial, 
respectively. 

In this work, we have made use of a recently proposed [TTl [IE] approach to 
quasiconformal mapping generation that does not need of functional minimizations. 
This method takes advantage of the fact that the transformation functions, x\x^ y) and 
y'{x^y)^ that minimize Winslow functional, also solve Laplace's equation: 

V^x' = 0; V'^y' = (18) 

where the Laplacian, taken in the physical space, is acting on the virtual coordinates. 
Therefore, solving Laplace's equation in physical domain subject to sliding boundary 
conditions, ft • Vx'\ leads to a quasiconformal map. Once the functions x'{x^y) and 
y'{x^y) are obtained, the gradients of the transformation can be computed in order 
to construct the inverse Jacobian matrix, = dx^' /dx\ By inverting it, the 

Jacobian matrix needed for the TO procedure is obtained. Now, assuming that the 
anisotropy of the quasiconformal mapping is negligible, we can approximate its metric 
tensor by a diagonal matrix and define the transformation medium by means of (15). 
Again, for problems where the polarization state of the EM wave is well defined an 
isotropic refractive index can be defined following (17). In this way, a quasiconformal 
transformation can be used to design an isotropic TO medium defined by the index 
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profile n = g^^^. 

3.1. Isotropic implementation of an SPP shifter 

As a first example of a plasmonic device based on isotropic dielectrics, let us consider the 
design of an SPP parallel shifter using a quasiconformal mapping. As in the previous 
examples, a TO device will be placed on top of the surface through which the SPP 
travels. In order to design an isotropic shifter with the same functionality as the 
anisotropic one previously considered, we start by defining the geometry shown in the 
inset panel of figure [6](b). The transformation medium, whose top view is shown in 
gray, is divided into three regions: an entrance for the SPP, of width //2, at the left, 
the shifting region (at 45°), of length rf, and an exit region at the right. Then, we solve 
Laplace's equations within the gray region with sliding boundary conditions ft • Vx' 
at the lower and upper boundaries and n • Vy' at the left and right boundaries. The 
resulting quasiconformal map is plotted in figure |6](b) , where the grid corresponds to 
lines of x'{x^y) = constant and y'{x^y) = constant. Following the method described 
above, we obtain an isotropic refractive index profile that varies between and 3.5 (color 
map in figure [6](b)). 




Figure 6. Simulation results of an isotropic SPP parallel shifter, (a) Transmittance 
through the shift for SPPs at different wavelengths. The red line shows 
the transmittance through a transformation medium characterized by anisotropic 
parameters (same as in figure |3|. The green and purple lines correspond to a parallel 
shifter characterized by the diagonal e^o and fiTO and the isotropic refractive index, 
respectively, both retrieved from the quasiconformal transformation. For all the cases, 
the TO parameters are implemented in the dielectric side only, (b) Quasiconformal 
map (grid) and refractive index profile (color map) of the device. The inset panel show 
the geometry layout, (c) Magnetic field, Hy^ of an SPP at 600 nm travelling through 
the isotropic shifter after implementing the index shown in panel b. 

The transmittance curve in figure [6](a) shows the results obtained from 3D 
simulations of the isotropic SPP shifter. SPPs at different wavelengths (from 800 to 1500 
nm) travel along a gold-vacuum interface and impinge on the transformation medium 
placed on top of the gold surface. As we have previously shown, a transformation 
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medium in the metal is not needed. The size of the device is the same as for 
the anisotropic shifter discussed above. Three different transmittance spectra have 
been calculated. The green and purple lines show the transmittance through the 
transformation medium resulting from the quasiconformal transformation. While the 



green line corresponds to the diagonal e^o and jiTo given by (15), the purple line 



corresponds to the index profile n{x^y) calculated according to (17) and shown in panel 
b. On the other hand, the red line shows the transmittance of the anisotropic shifter 
(same as in figurejs]), so that the efficiency of both approaches can be compared. Whereas 
the red line takes transmission values above 98%, the green and purple lines lie beneath, 
around 96% and 92% respectively. The are two main reasons for this feature. First, the 
fact that the transformation is quasiconformal leads to an impedance mismatch at the 
interfaces to free space, causing losses and lowering the green and purple curves. Second, 
the purple curve shows the lowest values in transmittance because when cto and jiTO 
are replaced by an isotropic index, the dispersion relation of the SPP is affected. In this 
case, the wave fronts of the SPP are distorted in order to follow the vertical lines of the 
quasiconformal map shown in panel b. This feature can be seen in figure [6](c), were an 
SPP at 600 nm is guided through an isotropic beam shifter. Despite the fact that an 
isotropic index reduces the transmittance, the shifting of the SPP is nicely accomplished. 

3.2. Isotropic implementation of an SPP right- angle bend 

A beam bend of the same characteristics as that discussed in section 2.3 can also be 
realized with isotropic materials. By means of a conformal transformation, a more 
realistic optical device that bends the SPP propagation direction by a right-angle can 
be designed. The anisotropic parameters (10), yielded by the transformation ([9]), are 



replaced by an isotropic refractive index if a conformal transformation is made [16j: 

n(r) = — (19) 
nr 

As illustrated in the inset of figure [7](a), the transformation medium, characterized by 
n(r), is a cylindrical shell, a < r < 6, of height hd placed on top of the metal surface. 
This refractive index smoothly decreases from 4.46 at r = a to 0.67 at r = 6. 

We have tested the functionality of the isotropic SPP bend by means of 3D 
simulations. The curves in figure [7|(a) compare the transmittance through the SPP 
bend with anisotropic parameters (red line) and bends characterized by parameters 
resulting from a conformal transformation (green and purple lines) over a wavelength 
range between 800 and 1500 nm. The green line shows the transmittance through a 
TO medium characterized by the diagonal cto and j^To calculated from the conformal 



transformation according to (15). The e^z and fi^z components of the tensors vary 



throughout the medium from to 16. The purple line corresponds to a bend 



characterized by the isotropic refractive index given by (19). In the last two situations 
the transmittance is limited to 96% due to impedance mismatch to free space. This 
feature, not present for the anisotropic case (red curve), arises when conformal 
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Figure 7. Simulation results of an SPP bend characterized by an isotropic refractive 
index. (a) Transmittance through the bend for SPPs at different wavelengths 
comparing bends with anisotropic, diagonal and isotropic parameters. The red line 
corresponds to the case of a transformation medium characterized by anisotropic 
parameters (same as in figure [4| . The green and purple lines stand for a conformal 
transformation: green curve for transmittance resulting from diagonal parameters and 
purple curve for the isotropic refractive index (top view plotted in the inset panel). 
In all cases, the TO parameters are implemented in the dielectric side only. The 
geometrical parameters are the same as in figure |4} (b) Electric field, E^, of an SPP 
at 1400 nm travelling through the isotropic bend. 



transformations are used. As shown in the inset panel of figure [7](a), the refractive index 
abruptly changes from n = 1 to n{r) at the interfaces along the entrance and exit faces 
of the bend. As a consequence, the SPP suffers refiections that lower the transmission. 
In this case, the green and purple lines lie very close to each other, meaning that the 
dispersion relation of the SPP is not affected by changing and jizz to n. The main 
reason for this is that the parameters resulting from the conformal transformation are 
smoother and lead to less distorted wavefronts than the quasiconformal transformation 
in the previous example of an SPP shifter does. Figure [7|(b) illustrates the performance 
of the bend for a gaussian SPP at 1400 nm. After impinging on the device, close to 96% 
of the SPP is bent in a right angle and leaves the transformation medium through the 
exit face to continue its propagation through the metal. 



3.3. Isotropic implementation of an SPP lens 

In this section we discuss the realization of an SPP lens with isotropic materials. The 
geometry and purpose of the device are the same as presented previously (see section 
2.4). An isotropic refractive index is obtained if, instead of the transformation given by 
( [TT] ), a conformal transformation between a semicircle and a rectangle is performed. In 
particular, we will make use of a Schwarz-Christoffel transformation as described in [16j. 
The refractive index profile that results from this transformation smoothly varies from 
to 1.31 and is shown in figure [8](b). To design the SPP lens a transformation medium 
with a rectangular cross-section characterized by the index distribution UToi^^y) 
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height hd is placed on top of a metal film. The device is surrounded by vacuum (n = 1). 
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Figure 8. Results obtained from 3D simulations on SPPs propagating through a lens, 
(a) Transmittance spectra through lenses characterized by different EM parameters. 
Red line: anisotropic cto and fiTO (same as in figure [5|. Green line: diagonal cto and 



ftTO resulting from a conformal transformation according to ( 15 ). Purple line: isotropic 
refractive index resulting from the same conformal transformation, (b) Refractive index 
profile yield by the conformal transformation, (c) Top-view of the y-component of the 
magnetic field. An SPP travels through a lens characterized by the refractive index 
shown in b and is focused within a width of the order of the wavelength. 

Figure [8] shows the results obtained from 3D simulations of SPPs propagating 
through the lens. The three transmittance curves in panel a compare the functionality 
of SPP lenses characterized by different EM parameters. The red line shows the best 
performance as it represents the transmittance through a lens characterized by the 
anisotropic and singular parameters given by ( [T2| ) (see figure [5]) . The green and purple 
lines correspond to lenses designed by means of a conformal transformation and therefore 
are situated at lower transmittance values due to impedance mismatch to free space. 
The green line accounts for the transmittance through a lens with diagonal parameters, 



cto and fiTO^ that are obtained from the conformal transformation according to (15). 
Only the third component of the tensors diagonal is different from 1, and it varies 
smoothly from to 1.72. The purple line corresponds to the isotropic refractive index 



retrieved from the conformal transformation according to (17) and is shown in panel b. 
The two cases considered for the conformal transformation (diagonal cto and j^To and 
isotropic n) give close transmittance results in this example, similarly to the SPP bend 
results above. Panel c shows the field pattern of an SPP travelling trough the isotropic 
lens as its wavefronts are focused on the exit face. 

3.4' SPP ground-plane cloak 

Finally, let us compare the two approaches to TO that we have presented (anisotropic 
and isotropic transformation media) as they apply to a ground-plane cloak for SPPs. 
This cloak was first proposed for mimicking a fiat refiective surface under plane wave 
illumination [15]. Given that TO is valid for any kind of EM waves, a ground-plane cloak 
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can also be used to suppress the scattering of an SPP from a bump on the metaUic surface 
along which the plasmon propagates (see figure [9]). Thus, the SPP appears to have 
travelled along a flat metal surface although it is physically travelling along the bump. 
Let us consider an SPP propagating in the x direction along a surface placed in the x-y 
plane. In order to cloak a bump deflned by the equation z{x) = ho cos^ (f^)? make 
a transformation in the x-z plane (perpendicular to the interface) from a rectangular 
region sitting on a metal surface to a region whose bottom boundary is lifted following 
the surface of the bump. We have considered two different transformations that yield 
the same functionality: a transflnite transformation, flgure|9](a), and a quasiconformal 
mapping, flgure[9](b). 

The transflnite transformation shown in panel a is made in the dielectric side of 
the interface in a box sized / x h and has an analytical expression (details can be found 
in [19j). Since it is not a conformal transformation, the EM parameters, cto and /xto? 
retrieved for the corresponding transformation medium are anisotropic. When the cloak 
is placed on top of the surface, an SPP smoothly travels across the bump without any 
scattering. Figure [9](c) renders the magnetic fleld component, Hy^ of an SPP at 700 nm 
travelling through a bump of height 0.2 //m. On the other hand, we have calculated 
the isotropic index proflle resulting from the quasiconformal map shown in panel b. In 
this case, the transformation medium is sized h' x V . The quasiconformal map was 
generated by solving Laplace's equation with sliding boundary conditions n-Vx' (at the 
lower and upper boundaries) and n • Vy' (at the left and right boundaries). The effect 
of a ground-plane cloak is achieved by setting = at the lower boundary so that 
the bump is mapped onto a flat surface. Then we have calculated the refractive index 
following the procedure described above. It ranges from 0.82 to L38 and the background 
value (n = 1) is recovered at the outer boundaries of the cloak. This means that the 
cloak is impedance-matched to vacuum, a fact that suppresses possible reflections at 
the interfaces. An SPP travelling through the isotropic cloak is illustrated in panel d, 
where the SPP smoothly passes the bump. In both cases (panels c and d), the wave 
fronts follow the transformed grid. In contrast, if no cloak is placed on top of the bump, 
the SPP suffers strong scattering as shown in flgure|9](e). 

4. Conclusions 

In conclusion, we have demonstrated that the application of the concept of 
Transformation Optics to Plasmonics brings into the fleld unprecedented control over the 
flow of SPPs. Based on this approach, we have presented the necessary tools to design 
feasible plasmonic devices. When considering SPPs propagating on metal-dielectric 
interfaces, a direct application of Transformation Optics would imply a manipulation 
of the optical properties of both the metal and dielectric materials. However, we have 
shown that modifying only the dielectric side leads to almost perfect functionalities 
because most of the SPP energy is carried within the dielectric medium at optical 
frequencies. Having proved this feature within a general framework in which the EM 
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Figure 9. Ground-plane cloak for a bump of height ho = 0.2 /im on a gold- 
vacuum interface infinite in the y direction, (a) Coordinate map of the transfinite 
transformation that leads to a ground-plane cloak characterized by anisotropic cto 
and fiTO' The cloak, of length / = 2 /im and height h = 2 /im, is placed in the 
dielectric side of the interface, (b) Refractive index profile (color map) resulting from 
the quasiconformal map shown by the grid. The size of cloak is given by = 6 /im and 
= 2.5 /im. (c) and (d) Magnetic field, Hy^ of an SPP at 700 nm travelling through 
the anisotropic (a) and the isotropic (b) cloaks, respectively, (e) SPP scattering on a 
gold bump when the cloak is not present. 



parameters are anisotropic, we have then presented plasmonic devices that can be made 
out of isotropic dielectric materials. Making use of conformal and quasiconformal 
transformations, control over SPPs flow can be achieved by placing a structured 
dielectric layer on top of the metal surface. Regarding practical implementations, this 
fact would render low loss and broadband performances. For instance, holes could be 
drilled in a dielectric layer to eflFectively mimic an isotropic refractive index. We expect 
that the flexibility of our methodology will ensure many diflFerent theoretical proposals 
and experimental realizations of plasmonic devices based on Transformation Optics. 
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Appendix A. Optical conformal mapping 



In the limit of ray optics, where the refractive index n{x^y) is nearly constant over the 
scale of a wavelength the field amplitude for each polarization satisfies Helmholtz's 

equation. In virtual space {z = x' + iy')^ this equation reads: 



V'2 + 



-n 



12 



= 



(A.l) 



where the Laplacian is V'^ = jdx'^ + jdx'^ . Let us now consider a transformation 
w = w{z) from virtual to physical space {w = x + iy) such that it is conformal. 
Making use of the relation between the Laplacian in virtual and physical spaces, 
we can transform Helmholtz's equation to physical space: 




dz 



duj 



n 



12 



7/; = 



This equation is left invariant if we choose a refractive index of the form [4J : 

dz 



n = 



duj 



n 



(A.2) 



(A.3) 



where n' is the background refractive index, no. By means of the Cauchy-Riemman 
equations we can express the refractive index in terms of the metric tensor of a conformal 
transformation: 

dx' 



n 



no 



d 
duj 



{x' + ly'] 



no 



dx dx 



no 



which is the same expression for the refractive index obtained according to TO (17). 
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